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We investigate aggregation driven by mass injection. In this stochastic process, mass is added 
with constant rate r and clusters merge at a constant total rate 1, so that both the total number 
of clusters and the total mass steadily grow with time. Analytic results are presented for the three 
classic aggregation rates Kij between clusters of size i and j. When Ki t j = const, the cluster size 
distribution decays exponentially. When Kij oc i + j or Ki,j oc i x j, there are two phases: (i) a 
condensate phase with a condensate containing a finite fraction of the mass in the system as well 
as finite clusters, and (ii) a cluster phase with finite clusters only. For Kij oc i + j, the cluster size 
distribution, Ck, has a power-law tail, Ck ~ k ' in either phase. The exponent is a non-monotonic 
function of the injection rate: 7 = r/(r — 1) in the condensate phase, r < 2, and 7 = r in the cluster 
phase, r > 2. 
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I. INTRODUCTION 

Aggregation processes in which small objects merge ir- 
reversibly to form larger clusters are ubiquitous in nature 
PI Q • For example, aggregation underlies the evolution 
of planetary systems in astrophysics Q , cloud formation 
and dust accumulation in atmospheric sciences 
as well as polymer and gel formation in chemical physics 
H, U 0. Aggregation also plays a central role in the 
theory of percolation , fractal formation , and net- 
work growth [l3l Il4| . 

Often, aggregation is driven by a constant injection 
of mass and consequently, the total mass grows indef- 
initely with time [ii GB [H E E HI- This is the 
case in chemical kinetics, and in particular polymeriza- 
tion where the aggregation rate depends on the number 
of clusters. Here, due to aggregation, the total number of 
clusters typically decreases with time or saturates at a fi- 
nite value. In such aggregation processes, the system may 
undergo gelation: a giant cluster develops and eventually, 
it contains all of the mass in the system 0, l22l |23| . 

In this study, we are interested in complementary ag- 
gregation processes that describe the growth of random 
structures such as random trees and rando m graphs, rele- 
vant in computer and information science |24ll25ll26ll27| . 
Here, the system consists of an ensemble of clusters where 
both the number of clusters and the total mass grow 
with time. For example, in the Internet, clusters are au- 
tonomous systems; injection models the creation of new 
autonomous systems and aggregation describes merger of 
different autonomous systems |2<Sl . 

We investigate aggregation processes with mass injec- 
tion where the total merger rate does not depend on the 
total number of clusters. Our main finding is that in 
such situations there is condensation rather than gela- 
tion. The system develops condensates that contain a 
finite fraction of the mass. These macroscopic conden- 
sates co-exist with microscopic clusters that contain the 



rest of the mass in the system. 

In our formulation, mass is injected at a constant rate 
and clusters merge at a constant total rate. We address 
the three classic kernels iCy for aggregation between 
clusters if size i and j, respectively: the constant rate 
Kij — const, the sum rate Kij =i+j, and the product 
rate Kij = i x j. 

For the constant aggregation rate, the cluster size dis- 
tribution decays exponentially with the cluster size. For 
both the sum and the product rates, where the aggre- 
gation rate grows with the aggregate size, the system 
undergoes a phase transition as a function of the injec- 
tion rate. When the injection rate is smaller than some 
critical value, the system is in a condensate phase. A fi- 
nite fraction of the total mass constitutes a macroscopic 
condensate, but the remaining fraction of mass is in the 
form of microscopic clusters. The condensate and the fi- 
nite clusters coexist. When the injection rate is larger 
than the critical rate, there are only finite clusters. 

In the cluster phase, the distribution of cluster size gen- 
erally decays as a power-law with the cluster size. How- 
ever, different behaviors emerge in the condensate phase. 
For the product rate, the size distribution falls-off expo- 
nentially at large size. For the sum rate, however, the 
behavior is always power-law. Interestingly, in the latter 
case, the decay exponent is not monotonic. It decreases 
monotonically with the injection rate in the condensate 
phase but increases monotonically in the cluster phase. 

The rest of this paper is organized as follows. We in- 
troduce the model, describe some of its basic features, 
and formulate the master equation approach in Sec. II. 
The constant rate, the sum rate, and the product rate 
are analyzed in sections III, IV, V, respectively. Gener- 
ally, our focus is the size of the condensate and the tail 
of the cluster size distribution. We conclude in Sec. VI. 
Technical derivations of two results in sections IV and V 
are presented in appendices lAl and iBl 
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II. THE MODEL 

In our model, there are two independent and com- 
peting processes: mass injection and merger of clusters. 
In the first process, monodisperse elemental clusters are 
added to the system. This injection process occurs at 
a constant rate. In the second process, two clusters are 
merged. The mass of the resulting cluster is equal to the 
sum of the two original masses. The total merger rate 
is constant as well, and since both processes occur with 
constant rates, we may set one of them to unity without 
loss of generality. We therefore set the mass injection 
rate to r and the merger rate to one. Also, since the in- 
jection is monodisperse, we set this injection size as the 
mass unit. Initially, the system contains no clusters. 

The total mass and the total number of clusters follow 
directly from these definitions. Each injection event in- 
creases the number of clusters by one and similarly, each 
merger event decreases the number of clusters by one. 
Thus, the average total number of clusters, N(t), satis- 
fies dN/dt = r — 1 and consequently, there is simple linear 
growth 



of the total mass, and summing the rate equations, we 
recover (J2J. 

We analyze the three classic aggregation rates: the 
constant aggregation rate Kij oc const, the sum aggre- 
gation rate Kij oc i + j, and the product aggregation 
rate Kij oc i X j 0, each representing a distinct grow- 
ing random structure. Generally, random structures such 
as random trees and random graphs are made of nodes 
interconnected by links. Mass injection represents ad- 
dition of isolated nodes, and merg er represents addition 
of a link between two nodes |25l |27| . In the constant 
rate process, two random structures are picked at ran- 
dom and an added link connects the two. In the product 
rate process, two nodes, picked at random, are connected 
by an added link. Finally, the sum rate is a hybrid of the 
constant and the product cases: the link connects a ran- 
domly selected node and a randomly selected structure. 
We note that the constant aggregation rate models an en- 
semble of random growing trees (no cycles are formed), 
while the sum and the product rates model an ensemble 
of random growing graphs. 



N(t) = (r-l)t. 



(1) 



We restrict our attention to situations where the number 
of clusters grows with time, r > 1. 

Merger events conserve the total mass, and hence, 
the mass changes only through injection events. With 
each injection event, unit-size mass is added to the sys- 
tem, and thus, the average total mass, M(t), obeys 
dM/dt = r. Consequently, the total mass also grows lin- 
early with time, 



M(t) =rt. 



(2) 



We investigate the cluster size distribution. Let Ck(t) 
be the average number of clusters of size k at time t. This 
quantity satisfies the master equation 



(3) 



i+j=k 



The initial condition is Cfe(O) = 0. The first term ac- 
counts for mass injection and the last two terms account 
for merger. The kernel Ki.j is defined as the aggregation 
rate between two clusters with size i and j. The total 
merger rate is constant, thereby implying the following 
constraint on the aggregation rate 



1 — 9 KijCiCj 



(4) 



III. CONSTANT AGGREGATION RATE 

For the constant aggregation rate, all pairs of clusters 
merge at the same rate, irrespective of their size. This is 
the simplest and the most widely used aggregation pro- 
cess with examples including fracta l aggre gates p^ . do- 
main growth [30| , and random trees |3lL l32| . We consider 
the following size-independent aggregation rate 
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(•5) 



This constant satisfies the normalization 0} . It decreases 
as the number of clusters increases so that the overall 
merger rate does not change with time. 

Substituting this constant rate into the master equa- 
tion ©, the cluster size density satisfies 



dC k 1 ^ 

i+j—k 



Ci Cj — Ck ■ 
N 



(6) 



The linear growth of the total density JTJ and the to- 
tal mass (0 suggest that the cluster size distribution 
also grows linearly with time. Indeed, the density of the 
smallest clusters obeys dC\/dt = r — 2C\/N so that this 
quantity, too, grows linearly with time, C\{t) ~ N(t)^j. 
Thus, we write the cluster size density as a product of 
the overall density N(t) and the time-independent cluster- 
size distribution 



The total density is of course N(t) — Ylk^k(t)- Sum- 
ming the master equation and using the constraint 
(QJ , we confirm the linear growth of the total density Q • 
Similarly, the total mass is M(t) = J^k kCk(t). Only the 
first term in the master equation affects the evolution 



C k {t)=N{t)c k . 



(7) 



The cluster size distribution is normalized, = 1. 

This form is consistent with the initial condition, and it 
satisfies © when the cluster size distribution obeys the 
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recursion equation 



with the constant 



(1 + r)c k = r <5 fei i + 



7 , °i C j ■ 

i+j=k 



(8) 



We utilize the generating function technique |33j to 
solve this equation. Let f(z) be the generating function 



fc=i 



(9) 



Normalization implies /(l) = 1. Multiplying Eq. © by 
z k and summing over k, the generating function satisfies 
the quadratic equation f 2 (z) — (l + r)f(z) + rz = 0. The 
solution is therefore 



/(*) = 



1 - 



4?- 



(1 



(10) 



We can confirm that f(z) = j^,z + • • • , in agreement 
with the above expression for C± . Writing the generating 
function as a power-series, the cluster size distribution is 
obtained explicitly, 



l + rr(fe-l/2) 
T(k + 1) 



4r 



(1 



(11) 



where T(x) is the Gamma function. Using the asymptotic 
property T(x + r)/T(x) ~ x r as x — > oo, we find 



Ck 

_ 1+r 



-3/2 



pk 



(12) 



with the constants a = -i±S= and 8 = Therefore, the 
cluster size distribution decays exponentially, but there is 
an algebraic correction. As the injection rate approaches 
the merger rate, r — ► 1, the cluster size distribution be- 
comes algebraic Ck ~ fc~ 3 / 2 . This limiting case coincides 
with the well-known behavior for time-independent con- 
stant aggregation rates [Tj| • 



IV. SUM AGGREGATION RATE 

Aggregation rates proportional to the sum of the clus- 
ter sizes are relevant in polymerization |34|, co agulation 
under shear- flows Q, and random graphs |2^. Subject 
to the constraint the sum aggregation rate is 



• J NM 

In this case, the master equation (0 becomes 

kN + M 



(13) 



dC k k m—^ 

= r dk i H > Ci C 

dt ' 2NM ^ 

i+j=k 



NM 



C k - (14) 



We again seek a solution of the form (JJJ. The cluster 
size distribution remains normalized, X)fc Cfc = 1' ano - ^ 
obeys the following recursion relation 



(k + R)c k = R5k,i + - Yj 



2 ^ ~ l 3 



(15) 



R = 



r - 1 



(16) 



This recursion relation can be manually solved to find 

ci = -j^i, C2 = ("R+i) ' e ^ c ' This procedure can 
be formally related to integer partitions following the so- 
lution procedure of Ref. |35|. Such a solution is useful 
only when the cluster size distribution decays sharply. 
Although it is difficult to obtain an explicit analytic so- 
lution as in (|ll(l . it is still possible to obtain many of 
the interesting properties of cluster size distribution from 
asymptotic analysis of the generating function. 

The generating function obeys the nonlinear ordi- 
nary differential equation 



*(/-!)£=*(/-.) 



(17) 



Derivatives of the generating function at z = 1 are re- 
lated to the moments of the cluster size distribution. For 
example, the average cluster size follows from the first 
derivative, (k) = /'(l). 

Differentiating (|17fl and then substituting z — \ yields 
a quadratic equation for the average cluster size 



R-^k) 2 - (k) + 1 = 0. 



(18) 



Naively, one expects that the average cluster size is the 
ratio between the total mass and the total number of clus- 
ters, (k) = ^ = ^-rj. Indeed, this is a solution of JTSJ. 
However, there is another solution (k) = r. This solution 
is not physical when r > 2 since the product N(k) can 
not exceed the total mass in the system. Each solution 
is relevant in the appropriate range of the parameter r, 
so that the full solution is 



(k) = 




(19) 



This assertion is supported by analysis below. The 
total mass contained by the clusters, M c = J2k^k, 
is of course proportional to the average cluster size, 
Mc = NY, k kc k = N (k)- Therefore, when r < 2, finite 
clusters contain only a fraction of the mass when r < 2, 
but they contain all of the mass when r > 2. 

This behavior can be reconciled with mass conserva- 
tion only if there is a condensate of mass M* that con- 
tains the remaining fraction = M*/M of the mass in 
the system as follows 



2 - r 1 <r < 2 
2 < r. 



(20) 



Thus the system undergoes a phase transition. When 
1 < r < 2, there is a condensate that contains a finite 
fraction of the mass. This condensate coexists with the 
finite clusters that contain the rest of the mass. The con- 
densate contains nearly all of the mass in the limit r — > 1: 
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FIG. 1: The exponent 7 versus r for the sum rate. 



when injection is very slow, the condensate "preys" on 
newly added mass. As the transition point is approached, 
the condensate mass vanishes, m* — > as r — > 2. When 
r > 2, the system contains only ordinary clusters. 

The tail of the cluster size distribution can be eval- 
uated from the z — > 1 behavior of the generating func- 
tion. In general, f(z) may contain both a regular compo- 
nent and a singular component, f(z) = f Ies (z) + /singO&)- 
The regular component is a power series in (z — 1). Let 
us assume a singular behavior with the leading behavior 



/ sing ex (1 - z)^ 1 



1, 



f( z ) = 1 + (k)( z - 1) + ■ ■ ■ + A{i - zy- 1 + 



(21) 



/reg(z) 



/stag (z) 



In the limit z — s 
tion Ijl7(l when 



1, this form satisfies the governing equa- 



<fc) 7 = R. 



(22) 



The algebraic form of the generating function implies an 
algebraic form for the tail of the cluster size distribution, 



Ck ~ k~ 



(23) 



as k — ► 00. The exponent is found by substituting H19|) 
into the relation H22fl . Therefore, there are two regimes 
of behavior 



7h Kr<2 
r 2 < r. 



(24) 



We note two remarkable features. First, the cluster size 
density is algebraic both in the condensate phase and in 
the cluster phase. Second, the characteristic exponent 
is a non-monotonic function of the injection rate: it de- 
creases monotonically with r in the condensate phase and 
it increases monotonically in the cluster phase (Fig. [I}. 

The exponent 7 is minimal, 7 = 2, at the phase tran- 
sition point, r = 2. For 7 < 2, mass conservation would 
be violated because the sum ^ fc kck is divergent. The 



restriction 7 > 2 justifies our previous choice of the non- 
trivial solution (k) = r in Eq. I|19fl . 

The behavior at the phase transition point requires a 
special treatment. We find that the cluster size density 
decays slightly faster than fc~ 2 , namely 



Ck 



k 2 (\nk) 2 



as 



(25) 



The derivation of this result is detailed in Appendix A. 
With this logarithmic correction, the sum ^ fc kck or the 
total mass in finite clusters, remain finite. 

In ordinary gelation (or percolation), the cluster size 
distribution is algebraic only at the critical point; away 
from criticality, it has exponential tails. In the present 
case, the cluster size distribution exhibits a strikingly 
different behavior — it is algebraic everywhere while at 
the critical point there is a logarithmic correction. 

The size of the condensate can be obtained directly by 
focusing on its dynamics. Let us assume that there is a 
giant cluster in the system with mass M 4 . Its growth rate 
is dM*/dt = Ki^i, i and substituting the aggregation 
rate lfH5|) . we arrive at 



dM* 
dt 



E 



N M 



M*(M-M*) 
JVM ' 



(26) 



In the second step, we made the approximation 
M* + i ~ M* as the condensate is much larger than the 
rest of the clusters. When the condensate contains a fi- 
nite fraction of the mass, M* = m*M, then from 12611 . 
we find m*(r — 1) = m*(l — m*). Solving this equation, 
we recover the condensate mass (|20fl . 

This approach can also be used to derive the size of 
the largest cluster in the cluster phase. From (|20|l . the 
quantity M* is negligible compared with M, and there- 
fore, the governing equation (|26|l reduces to dM*/dt w 



(r-l)t 



M so that 



(27) 



with 7 = r when r > 2. Therefore, the size of the largest 
cluster grows algebraically with time. This result can be 
alternatively derived using the algebraic behavior IL'MI) 
and the extreme statistics criterion 1 ~ N^2 k>M c&. 

V. PRODUCT AGGREGATION RATE 

The product aggregation rate models polymerization 
and gelation H, U as well as random graphs 0> E3- 
and in our case it has the following form 



1X1 

Ka = 2 — 4- ■ 



(28) 



The explicit rate equation for the product aggregation 
rate is 



(29) 



i-\-j—k 



Since the total mass appears in the denominator in the 
rate equation, we use a different normalization 
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C fe = Mc k . 



(30) 



The average number of clusters of size k still grows lin- 
early with time. The transformation 1)30(1 reduces the 
master equation ((29(1 to the nonlinear recursion equation 



rcfc = rS k i + 



E 

i+j—k 



ij Cj Cj - 2kc k . 



(31) 



Given of the structure of this equation, it is convenient 
to use a different definition of the generating function 



/(*) 



k=l 



kckZ 



(32) 



The generating equation satisfies the very same equation 
lfT7)l with 



R 



(33) 



Although the governing equation is the same, the bound- 
ary condition is different. Since the product aggregation 
rate grows with the cluster size, we expect that again 
there are two phases: a condensate phase and a clus- 
ter phase. The total mass contained in finite clusters 
is given by M c = M/(l). Therefore, in the condensate 
phase /(l) < 1 while in the cluster phase /(l) = 1. This 
change in the boundary condition results in a drastically 
different behavior. 

We first discuss the cluster phase where /(l) = 1, and 
consequently, the analysis is a straightforward general- 
ization of the above. The first derivative at z = 1 again 
satisfies i? _1 [/'(l)] 2 - f'(l) + 1 = 0. At the critical point 
R c = 4 and therefore r c — 8. Solving this quadratic 
equation, the first derivative is 



/'CD = i 1 



r 



(34) 



The first derivative is now the ratio between the second 
and the first moments of the cluster size distribution, 
/'(l) = (k 2 )/(k). The tail behavior follows from the 
singular component of the generating function as in 121|) 



/(*) - 1 + f(i)(z - 1) + • • • + b(i - z y- 2 + 



(35) 



/=ing(z) 



This again implies the power-law decay (|23|l for the clus- 
ter size distribution. We note that the unit shift in the 
exponent is due to the different definition of the gener- 
ating function lj3*2"l) . From the governing equation (|T7jl . 
the exponent 7 satisfies (7 — 1)/'(1) = R, and the decay 
exponent is 



7 = 1 



(36) 
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FIG. 2: The exponent 7 versus r for the product rate. 



Similar to ((24(1 , the characteristic exponent grows linearly 
with the injection rate, 7 ~ |, at large injection rates 
(Figure [21). The minimum value, 7 = 3, is achieved at 
the phase transition point. A more careful analysis is 
again required at the phase transition point; we find (the 
analysis is essentially the same as in the case of the sum 
rate, Appendix lA"|l that c k ~ k~ 3 (\nk)~ 2 as k — > 00. 

In the condensate phase, we are able to obtain the 
mass of the condensate only in the vicinity of the phase 
transition, 



exp 



tta/8 



(37) 



as r i 8. The derivation of this result resembles that of 
Refs. [HE^; it is detailed in Appendix FBI In contrast 
with the sum rate, the phase transition is now very gentle: 
all derivatives of the condensate mass vanish at the tran- 
sition point r = r c 38, 39]. In practice, it may be difficult 
to locate the phase transition point (Figure 0). We com- 
ment that similar behavior was rece ntly found in several 
models of growing networks [M 03 Sfl El El El E| 

We performed Monte Carlo simulations to obtain the 
condensate mass shown in Figure [3] In the simulations 
there are two elemental steps: injection with probability 
r/(r + 1) and aggregation with probability l/(r + 1). In 
an injection step, a cluster with unit mass is added into 
the system. In an aggregation step, two clusters, picked 
with probability proportional to the product of the two 
masses, are merged. The simulations results represent a 
single run with a total mass of M = 10 6 . 

As a signature of the phase transition, the quantity 
/'(l) has a discontinuity at the phase transition point. In 
the condensate phase, this quantity is obtained, using the 
fact that /(l) < 1, directly from the governing equation 
(|T7jl . /'(I) = R = r/2. In the cluster phase it is given 
by 134(1 . Therefore, the ratio between the second and the 
first moments has a jump at r c = 8, 



(k 2 



(k) 



r I r c 



(38) 
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FIG. 3: The condensate mass m* as a function of r obtained 
from Monte Carlo simulations of the aggregation process. 



A similar jump occurs in the Berezinskii-Kosterlitz- 
Thouless phase transition [3j| ^ij^ . 

Finally, we study the tail of the cluster size distribu- 
tion. For this purpose, it is convenient to make the trans- 
formation z — e w so that the generating function i|32|) is 
redefined 



E 

k=l 



k Cfce 



kw 



(39) 



Substituting this definition into the recursion equation 
(|31[) . the generating function obeys 



df _ f-e w 
dw / - 1 



(40) 



The denominator suggests that f(w) has a singularity at 
wo with f(wo) = 1. From the definition (|39[1 . this implies 
the exponential decay Ck ~ exp(— k/ko) with ko — 1/tWo- 
Since we can not solve for the generating function, wc 
can not locate this singularity explicitly. Nevertheless, 
one can still deduce the behavior near this singularity 
using asymptotic analysis. Again, we assume that the 
generating function has a regular component and a sin- 
gular component near w = wq, 



fH 



i+_^+ A{w - w y~ 2 

/rcgO) 



(41) 



Substituting this form into (|40|l . and equating powers of 
(wq — w) on both sides of the equation we obtain v = 5/2. 
Using the definition l|39fl . the leading behavior of the 
singular component 1 — f(wo) ~ (wq — w) v ~ 2 implies 
an algebraic correction to the leading exponential decay, 
Ck ~ k~ u exp(— k/ko) as k — > oo. We conclude that the 
tail of the cluster size distribution decays as follows 



Ck 



(42) 



The characteristic size can be related to the conden- 
sate mass in the vicinity of the phase transition as the 



characteristic size is expected to diverge in this limit, 
ko — * oo. Consequently, the singularity is located close 
to the origin, wq — > 0. Let us estimate the behavior of 
the generating function near the origin. From the defini- 
tion 139fl . we have f(w = 0) = 1 — to*. Also, from the 
governing equation we have /'(0) = 4 in the limit r f 8 
as in J3HJ. Therefore, f(w) = 1 — m* + 4w as w — > 
and from the condition f(wa) — 1, the location of the 
singularity is wq = m*/4. Therefore, the characteristic 
size grows as follows 



feo ~ 1/m* 



(43) 



in the vicinity of the phase transition point, r f 8. 

The tail behavior coincides with the critical behavior 
at sufficiently small sizes, Ck ~ k~ 3 for k <C fco, and 
exponential decay, Ck ~ exp(— k/ko), occurs beyond that 
scale, k^> ko. The two behaviors should of course match 
at k w ko'- this implies the proportionality constant in 



63, 



Ck 



-1/2 



k~ 5 / 2 exp(-fc/fc ) 



VI. CONCLUSIONS 



In conclusion, we studied aggregation with constant 
injection of mass. In this process, the total number of 
clusters grows with time. For aggregation rates growing 
as the sum or the product of the cluster sizes, there are 
two phases: a condensate phase and a cluster phase. In 
the condensate phase, a condensate containing a finite 
fraction of the mass coexists with finite clusters, while in 
the cluster phase there are only finite clusters. 

For the sum rate, the mass of the condensate is a linear 
function of the injection rate. Also, the cluster size dis- 
tribution decays algebraically in both phases and inter- 
estingly, the decay exponent is a non-monotonic function 
of the injection rate. For the product rate the conden- 
sate mass is extremely small in the vicinity of the phase 
transition point and consequently, the phase transition 
is very gentle. In this case, the tail of the cluster size 
distribution is exponential in the condensate phase but 
algebraic in the cluster phase. 

We comment that there are two frameworks for de- 
scribing aggregation processes: the Flory approach that 
allows interaction between giant and finite clusters @ 
and the Stockmayer approach that allows for interactions 
between finite clusters only 0. We used the more chal- 
lenging former approach as it is the appro pria te approach 
for modeling growing random structures [45| . 

The various aggregation processes correspond to dif- 
ferent random growing structures, but our study focused 
only on the size of these structures. We note that this 
theoretical framework can be generalized to also study 
structural properties such as paths and cycles |27J. 

There are a number of possible extensions of this work. 
We focused on the three classic aggregation rates where 
the generating function obeys closed equations. This 
framework does not allow derivation of the necessary con- 
ditions for the emergence of a condensate as a function 
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of the aggregation rate. Based on the sensitive algebraic 
behavior in both of the phases, we speculate that the sum 
rate may be the marginal case for condensation. 
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APPENDIX A: DERIVATION OF J25l) 

At the phase transition point r = 2 we have 
7 = 2 and therefore the leading singular term 
/sing oc (1 — z) 7-1 becomes regular. This suggests to use 
instead /sing oc (1 — z)u(z) where u(z) vanishes slower 
than any power of (1 — z) as z — > 1. Thus at the phase 
transition point we employ the following expansion 

f(z) = 1 + 2(2-1) + (z-l)u{z) + ... (Al) 

Substituting QAlj l into l|17|) yields the differential equa- 
tion 

„ du u 2 , . , 

( z - 1 T + V— = ( A2 ) 
dz 2 + u 

whose (implicit) solution is 

2 

hlnu + khl - z) = const. (A3) 

u 

In the limit z — ► 1, the integration constant is negligible 
compared with the logarithmic term in the limit z — > 1 



8 



and consequently, u — > 2/ln(l — z). Indeed, u vanishes 
slower than any power of (1 — z) as z — > 1. Thus 



f(z) = 1 + 2(^-1) + 



2{z-l) 



(A4) 



ln(l-z) 

Inverting this expansion leads to Eq. 36] . 

APPENDIX B: DERIVATION OF 11571) 

The mass of the condensate follows from the behavior 
of the generating function at z = 1, m* = 1-/(1). 
To analyze the behavior near this region, we make the 
transformations 



f(z) = l + xg(x), 
z = 1 — x. 



(Bla) 
(Bib) 



With these transformations, the equation for the gen- 
crating function l|17|) is transformed into the following 
first-order nonlinear differential equation 



x gg' + g 2 + Rg + R = 0. 



(B2) 



In writing this equation, we kept only the leading order 
terms. Writing g2+ ^ g+fi c?ff + ~dx=0, and integrating, 
we have 



- ln(g + Rg + R) + In x tan 

2 2a 



const 



where a = \/R — R 2 /4. The integration constant 
can be evaluated by taking the x — » limit. Using 
to* = — lirnx_>o xg{x), the first two terms in the above 
equation approach In to, in the limit x — > 0. Using 
lirn^o g(x) = — oo, the last term approaches Hence 



1 



ln(# 2 + Rg + R) + In a; 



2a V a 



In i 



4a " 



(B3) 



Since we are interested in the behavior near the phase 
transition point, we take the limit r — > 8. In this limit, 
we can replace R by 4 and also, the quantity g 2 + Rg + i? 
by (<? + 2) 2 . Additionally, we may replace tan -1 ^i 2 - by 
§ — -^2. With these substitutions, Eq. I)B3|I becomes 



ln(-g - 2) + In: 



2 27T 

= In to, H . (B4) 

g + 2 a 



Next, we evaluate the left-hand side precisely at the 
phase transition point, r = 8. The critical behavior is 
detailed in Appendix A. Substituting g = —2 — u and 
x = 1 — z into HA3fl . gives 



ln(-g - 2) + h\x - 



const. 



g- 



(B5) 



Substituting this into (|B4|I we obtain the condensate 
mass in the vicinity of the phase transition l|37() . 



